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Abstract 

The analysis of the grey tracks produced in Deep Inelastic Scattering (DIS) off nu- 
clei provides important information on the space-time development of hadronization 
in nuclear medium. This method is complementary to the measurement of nuclear 
attenuation of leading inclusive hadrons. While the latter is focused on the hadroniza- 
tion dynamics for the quite rare process of leading hadrons production, the former 
covers the main bulk of inelastic events, and its Q 2 dependence is a very sensitive 
tool to discriminate between different models of hadronization. Employing the model 
of perturbative hadronization developed earlier, we calculate the Q 2 and XBj depen- 
dences of the number of collisions and relate it to the mean number of grey tracks, 
using an empirical relation obtained from the analysis of data from the Fermilab E665 
experiment on DIS of muons off the Xe nucleus. We found the number of grey tracks 
to rise steeply with Q 2 in good agreement with the experimental data. 
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Nuclear targets serve as a natural and unique analyzer of the space-time development 
of strong interactions at high energies PQ |2J. Due to Lorentz time dilation, the projectile 
partons may keep their coherence for some time, but once they become incoherent, the 
cross section of the final state interaction (FSI) increases and the nucleus may act as a 
filter of hadronization mechanisms, provided quantities which are sensitive to the nuclear 
modifications of FSI are measured. 

One way for experimental testing of various theoretical models which have been proposed 
to investigate the space-time evolution of hadronization, is the measurement of the nuclear 
modification factor in inclusive production of leading hadrons [HI EJ |3] • Recent experimental 
data from the HERMES experiment at HERA [5] have confirmed the predictions made in 
PP, based on the perturbative gluon radiation approach. Attenuation of the leading hadron 
is partially ascribed to the FSI of the pre-hadron (a colorless qq dipole prior formation of the 
hadronic wave function) produced in the color neutralization of the radiating quark which 
was created in the hard 7*-nucleon interaction. However, other approaches for in-medium 
hadronization have recently been employed, which are also able to provide a reasonable 
explanation of the HERMES data. We mean, in particular, the approach of Ref. jB], in 
which nuclear attenuation is explained in terms of the induced energy loss scenario, in 
which the hit parton loses energy by hard multiple scattering with other nucleon's partons, 
and the approaches of Ref. [3 El El, based upon the color string model. Thus it appears that 
other, more selective phenomena which could discriminate various models of hadronization 
are required. In this respect, the investigation of the Q 2 dependence of nuclear effects would 
probably be one of the best candidates 1 . As a matter of fact, recent data from HERMES 
are in a good agreement with calculations performed in [2] , but in strict contradiction with 
the Q 2 dependence predicted in [SJ. Unfortunately the HERMES data exhibit only a mild 
Q 2 dependence and therefore could not be very selective for other hadronization models. 

Note that production of leading hadrons with large fraction Zh of the initial parton 
momentum is a rare, nearly exclusive process, which has quite a specific time development. 
In the main bulk of events, the jet energy is shared by many hadrons and therefore in this 
case it is a difficult task to find observables which are sensitive to the time development of 
hadronization. In Ref. [TDj it has been proposed as a possible candidate the DIS exclusive 
process A(e, e'B)X, where B is a detected recoil nucleus with a less number of nucleons 
than the target one [TT], and X, as usual, denotes the unobserved jets of hadrons resulting 
from hadronization. It has indeed been shown JU| that the survival probability of the 
recoil nucleus B in such a process is rather sensitive to the time evolution of hadronization 
and, particularly, to its Q 2 dependence. This quantity, which is quite difficult to access 
experimentally in case of heavy nuclei, could be investigate using few-body nuclei. Indeed, 
first experimental data were recently obtained at Jefferson Lab ^2] for a deuteron target, 
i.e. by measuring the recoiling protons in the DIS process D(e,e'p)X in a wide range 
of kinematics. Preliminary results of calculations J3| show a good agreement with the 
predictions of the hadronization model of Ref. p]] . 

It should however be pointed out that the dominant channels of DIS are associated with 
events where the recoil nucleus B does not survive, but breaks apart to fragments, whose 

this note Q 2 = —q 2 = q 2 — v 2 denotes the squared four-momentum transfer 
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distribution in Q 2 and other kinematic variables could also shed light on the time evolution of 
the jet in the nuclear matter. A typical process of this type is the production of the so called 
"grey tracks" (GT) in DIS semi-inclusive processes. GT are those hadrons, predominantly 
protons, whose momenta are in a few hundred MeV/c range. They have been observed in 
various DIS processes induced by different projectiles; in this note we will consider the GT 
production in the Fermilab E665 experiment on [i — Xe and fi — D processes at 490 GeV 
beam energy; in this experiment the GT have been associated to protons with momentum 
in the interval 200 - 600 MeV/c. 

The investigation of the Q 2 -dependence of the GT production may bring forth precious 
information about multiple FSI of the jet originated from the DIS in the nuclear medium. 
The mean number of GT, to be denoted henceforth as (n gr ), correlates with the intensity 
of the FSI: this is why the production of GT has been traditionally used as a probe for the 
dynamics of the interaction and the centrality of collisions. The latter is usually character- 
ized by the mean number of collisions, (z/ c ), i.e. by the number of bound nucleons which 
took part in the interaction 2 . 

The dependence of (n gr ) upon Q 2 , v and xsj has been measured in the E665 experiment 
and it is the aim of this note to present a theoretical interpretation of such a dependence 
resulting from a specific model of hadronization. To this end we will proceed in the following 
way: our hadronization model allows one to calculate the average number of collisions (u c ) 
which is related to (n gr ). The relation between the observed (n gr ) and models for (u c ) 
can be obtained with various Monte-Carlo generators for hadron cascading in nuclei. The 
results of these calculations however do not seem to be very accurate and a more reliable 
and model independent approach has been adopted in , where the following relation has 
been found, 

(u c ) = (2.08 ± 0.13) + (3.72 ± 0.14) (n gr ) , (1) 

basing on the measured average total hadronic net charge and its relation to (u c ). 

We have used the same relation Eq. (JTJ) to obtain (n gr ), after calculating the average 
number of collisions (u c ) occurring during the evolution of a jet produced in DIS and prop- 
agating through nuclear matter. For minimum bias DIS events we obtain, 

oo oo 

(u c ) = J d 2 b J dzp A (b,z) J dz' p A (b,z')a eff (z-z') + 1 . (2) 

— oo z 

This equation results from the following description of the DIS process: the hard interaction 
of the lepton takes place on a bound nucleon at impact parameter b and longitudinal 
coordinate z; then the knocked out quark hadronizes and initiates a jet which propagates 
through the nucleus interacting with other bound nucleons with the effective cross section 
a e ff(z — z') depending on the distance (or time, provided that the quark propagates with the 
speed of light). The first term in Eq. (j2J) gives the amount of new nucleons involved in the 
process, i.e. the number of collisions, whereas the second term represents the contribution 
from the recoil nucleon formed in the initial, hard collision between 7* and one of the bound 
nucleons. 

2 Note that this number is different from the expansion parameter, <Ji n TA(b) in the Glauber formula 
which is also frequently called number of collisions (see the discussion in 
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Evaluation of Eq. (J2J) requires an explicit form for the time-dependent effective cross 
section <J e ff, which can only be obtained within a model for hadronization. To this end, 
we employ the model suggested in JU] : it combines the soft part of the hadronization 
dynamics, described in terms of the string model, with the hard part, described within 
perturbative QCD. There are many experimental evidences showing that gluons are located 
within small clouds around the valence quarks fTf\ ITH] . Therefore, if Q is less than the 
mass scale A = 0.65 GeV controlling the transverse quark-gluon separation, they can hardly 
be shaken off; in this case the string model is a proper description of hadronization. If, 
however, Q > A, perturbative gluon radiation should be taken into account. Employing 
A as the bottom limit for the integration over the gluon transverse momenta (see below), 
double counting is excluded when the string and pQCD contributions are added to describe 
hadron production. This model of hadronization is close to the one used in PjJ |2] , where the 
string contribution was mocked up by the soft part of gluon spectrum, and the predictions 
of the two versions of the model are rather similar. 

It was found in ^U] that the multiplicity of the produced hadrons, or better to say the 
pre-hadrons which are colorless dipoles, rises with time as 

(n h (t)) =n M (t)+n G (t) , (3) 

where n M is the amount of the pre-hadrons produced due to decays of the string, and na{t) 
the one produced by gluon radiation. 

The string contribution to (rih{t)) has been found in employing the standard dy- 
namics of string decay OI^D]- The pre-hadron multiplicity as function of time, n,M(t), was 
found in the following form 10 , 

, , In(l + t/A£) 

Here the time scale At is related to the probability w of a qq pair creation in the color 
field of the string, per unit of time and per unit of string length. This parameter evaluated 
either in the Schwinger model or from the phenomenology turns out to be w ~ 2fm -2 . 



Correspondingly, At = y2/w = lfm. 

Note that the logarithmic dependence on time in Eq. (jlj) is rather obvious. It is related 
to the fact that hadrons produced via string decays build a plateau in rapidity. Since the 
momenta acquired by pre-hadrons are proportional to the time taken by the string decay, 
the number of decays rises linear in the ln(t) scale. 

Concerning the gluon contribution to the hadron multiplicity Eq. (|3*|). we employ the 
large N c approximation and replace each radiated gluon by a color octet qq pair, and then 
combine the quarks and antiquarks into colorless dipoles (pre-hadrons). This is the origin 
of the second term na(t) in (J3J). 

It is well known that radiation does not happen instantaneously, but takes time called 
coherence time, t c = 2oo/k^, where u and kx are the energy and transverse momentum of 
the radiated gluon. At shorter times the quark-gluon system is still in coherence with the 
initial quark. Integrating the perturbative gluon radiation spectrum over u and k^ with a 
weight factor Q(t — t c ) one gets the amount of gluons which lost coherence, i.e. were radiated 
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over time interval t after the DIS interaction. The time dependence of the gluon radiation 
was found in [TU] to be controlled by the parameter £ = (m^ %Bj)~ l = 0.2im/xsj, and the 
number of perturbatively radiated gluons was evaluated in the following form 



n G (t) 




+ In 



tA 



qcd \ 



In 



HQ/AqCd) 

ln(A/ Aqcd 



(5) 



for t < t . At longer times, t > t , the t-dependence starts leveling off, 



n G {t) 



+ 




+ In 



tA 



f 2 ) ln 



HQ/A Q cDyJto/t) 

\n(\/A QCD ) 



HQ/A QCD ) 



\n(Q/ A Q cd Jt /t) 



(6) 



and saturates at t > t Q 2 /X 2 = 2z//A 2 , which is a very long time interval, substantially 
exceeding the nuclear size for the energies v = 50 — 400 GeV covered in the E665 experiment. 

Having obtained the explicit forms of the string and perturbative radiation contributions 
in Eq. (J3J), let us now discuss the interaction cross section of the produced colorless dipoles. 
These, as we mentioned, should be treated as pre-hadrons, until they form hadronic wave 
functions and get definite masses. As for their interaction cross section, what matters is 
their transverse size, rather than their mass. Since the decay of a string is a soft process, it 
is natural to assume that the decay products have a size of the order of the hadronic size 
and therefore interact with about the same cross section. 

Less obvious is the situation with pre-hadrons produced from perturbative gluons. Since 
the mean transverse momenta of the gluons follow the photon virtuality Q 2 , the initial size 
of the produced pre-hadrons should be small, of the order oil/Q. Then one may expect the 
color transparency effect to be at work j2U E2] • Indeed, this happens for leading hadrons, 
and the evolution of the pre-hadron size has a strong impact on nuclear transparency PQEj- 
In the present case, however, we deal with the slowest part of the hadron spectrum which 
is independent of the photon energy u, provided that it is sufficiently high. Only those 
pre-hadrons are produced within the nuclear volume whose energy does not exceed Eh < 
Ra (Q 2 — A 2 )/2 ln(Q 2 / A 2 ). The size of such a pre-hadron evolves so fast that the largest part 
of the effect of color transparency is washed out. Even for the most energetic exclusively 
produced hadrons in the E665 experiment, the effect of color transparency is rather moderate 
it is much weaker and hardly detectable at lower energies [25J. Therefore, we neglect 
these corrections and assume that pre-hadrons interact with the hadronic cross sections. 

Eventually, we are in the position to evaluate the time dependent effective absorption 
cross section in (J2J), to be used in the description of GT production; by treating all the qq 
colorless dipoles as mesons (M), we obtain 



eff 



(t) 



MN 



n M {t) +n G (t) 



(7) 



Note that this expression is different from the one used in [10]. Firstly, we excluded the 



term corresponding to interaction of the nucleon originated from the first DIS event, 
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since we treat this nucleoli as a participant giving rise to the last term in (J2J). Secondly, we 
replaced the total meson-nucleon (MN) cross section by the inelastic one. This is because 
most of recoiling protons in elastic and diffractive scattering have too small momenta to 
contribute to the production of GT. Therefore, we use a^ N = — a^ N — = 17.7 mb. 

Before illustrating the calculation of the Q 2 dependence of the mean number of GT 
n gr(Q 2 ), let us briefly discuss the expected dependence of n gr upon the energy transfer v. 
Both the string model and perturbative radiation lead to the mean multiplicity which rises 
logarithmically with energy. Naively, one may conclude that the amount of GT follows the 
multiplicity and should rise with energy as well. However, as already discussed, it takes 
time to produce hadrons. The observed multiplicity does not emerge in an instantaneous 
explosion-like particle production, but is a result of time development of the hadronization 
process. The full time of jet production is proportional to its energy and exceeds the nuclear 
size at the E665 energies. However, only those gluons which are radiated inside the nucleus 
can contribute to the production of GT. With the condition t c < one should integrate 
the radiation spectrum from u min = k T up to u) max = k\R&- The result is independent of 
the jet energy. Indeed, this is confirmed by the results of the E665 experiment (see Figs. 2, 
17 in [TU). 

At the same time, one should expect a Q 2 dependence of the mean number of GT. 
Indeed, the stronger kick gets the quark, the more gluons it shakes off. The mean transverse 
momentum of radiated gluons follows Q 2 which plays the role of the upper cut-off for 
the integration (Pkr/k 2 ^. Correspondingly, the multiplicity of the radiated gluons rises 
logarithmically with Q 2 . 

Let us now quantitatively analyze the behavior of (n gr (Q 2 ))- To this end we have calcu- 
lated the mean number of collisions by using Eqs. (J2J) and (J7J). Then, applying (0), we have 
obtained the number of GT, (n gr (Q 2 )). We fixed XBj at the mean value (xzj) = 0.068 [Tlj . 
The result shown by dashed curve is compared with the experimental data from the E665 
experiment (Hj in Fig. ^ 

If the time interval to is of the order or smaller than the nuclear size, the chosen value 
of xbj affects hadronization and is important. The E665 data are subject to quite a strong 
correlation between Q 2 and XBj, which may change the results depicted in Fig. [T] We 
parametrize this correlation as Q 2 = A + B XBj, where the parameters A = 2.2 GeV 2 and 
B = 178 GeV 2 are found using the values of (Q 2 ) and (xsj) for two regions called in [Tl] 
"shadowing" and "non-shadowing". At very small XBj < 0.01 we introduce an additional 
XBj dependence of parameter A which is, however, unimportant for the kinematic region 
under discussion. Our results corrected for the Q 2 — XBj correlation are depicted by the 
solid curve in Fig. [T] In spite of the difference between the two curves, both agree with the 
data. 

In order to see the dependence of (n gr ) on xsj explicitly we calculated it either fixing Q 2 
at the mean value (Q 2 ) = 14.3 GeV 2 , or applying the Q 2 — XBj correlation introduced above. 
The results are depicted in Fig. by the dashed and solid curves, respectively. Although 
both curves slightly overestimate the data, the solid one better reproduces the measured 
slope. 

In closing, the following comments are in order: 

• In our calculations we employed the same hadronization model as in QjJ |2J EI] with 
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Figure 1: The mean number of grey tracks (n gr ) produced in the [iXe 
DIS as function of Q 2 . The theoretical predictions based upon Eqs. (JTJ), 
(0) and (J7J), are compared with the experimental data from Ref. [TI] . 
which correspond to non-shadowing (xBj > 0.02) region. The dashed 
curve corresponds to fixed XBj = (xBj) = 0.068. The solid curve includes 
the correlation between Q 2 and XBj (see text). 



no readjustment of the parameters, and the empirical relation between the number 
of grey tracks and the number of collisions given by Eq.Q. The model successfully 
passed the new test; 

• the results of calculations may be subject to further corrections. In particular, we 
neglected resonance decays, which is justified by their rather high energies. On the 
other hand, the formation time of hadronic wave function is much shorter than the 
resonance life time and is rather short for pre-hadrons produced within the nucleus. 
This is the reason why we neglected the effects of color transparency which are impor- 
tant for leading hadrons [HI2|- This may explain why our results exhibited in Fig. ^ 
somewhat overestimate the data; 

• the Q 2 dependence of the mean number of grey tracks serves as a sensitive tool for test- 
ing theoretical models of hadronization. It can be seen that the amount of grey tracks 
doubles within the range of Q 2 covered by the kinematics of the E665 experiment, and 
our calculations well reproduce such a steep variation of (n gr (Q 2 ))] 

• whereas nuclear attenuation of leading hadrons carries information of space-time devel- 
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Figure 2: The mean number of grey tracks (n gr ) produced in the [iXe 
DIS as function of Xsj- The experimental data are from Ref. [HJ. The 
dashed and solid curves show the results of calculations with fixed Q 2 = 
14.3 GeV 2 and corrected for the Q 2 — xgj correlation respectively. 



opment of hadronization in rather rare events, grey tracks provide precious information 
about hadronization dynamics for the main bulk of DIS events; 

• We found a rather flat dependence of (n gr ) on xsj at fixed Q 2 . However, the Q 2 — XBj 
correlation existing in the E665 experiment leads to a growth of the number of GT 
with xsj- Our calculations corrected for this effect well describe the measured XBj 
dependence of GT; 

• in our calculations we assumed that the value of Bjorken scaling variable xsj is suf- 
ficiently large to neglect the sea, and this is the reason why we compared our results 
with the E665 data at large values of x B j. At small values of xbj one should take 
into account the production of two jets with different momenta, with only one of them 
producing a grey track corresponding to the target nucleon on which the DIS process 
occurred. On the other hand, the E665 experiment detected no substantial change 
in behavior of different observables in the shadowed compared to the nonshadowed 
regions. This was explained by a large contribution from diffraction [Tlj . 
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